Dynamic effective mass of granular media 
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We develop the concept of frequency dependent effective mass, M(uo), of jammed granular 
materials which occupy a rigid cavity to a filling fraction of 48%, the remaining volume being 
air of normal room condition or controlled humidity. The dominant features of M{uj) provide 
signatures of the dissipation of acoustic modes, elasticity and aging effects in the granular 
medium. We perform humidity controlled experiments and interpret the data in terms of a 
continuum model and a "trap" model of thermally activated capillary bridges at the contact 
points. The results suggest that attenuation in the granular materials studied here can be in- 
fluenced significantly by the kinetics of capillary condensation between the asperities at the contacts. 

PACS: 45.70.-n, 46.40.-f, 81.05.Rm 



A distinct feature of jammed or loosely packed gran- 
ular materials made of a variety of different materials 
such as sand, steel, polymer or glass is the ability to dis- 
sipate acoustic energy through the network of interparti- 
cle contacts or viscous dissipation through the surround- 
ing medium. Indeed,loose grains damp acoustic modes 
very efficiently 0, 0, H, HI, and they are routinely used 
as an effective method to optimize the damping of un- 
wanted structure-borne acoustic signals [l[. Despite its 
fundamental importance and practical applications, the 
microscopic origins of the mechanisms of dissipation in 
jammed granular materials are still unknown. 

In this Letter, we pursue the concept of a frequency 
dependent effective mass, M(u>), of a loose granular ag- 
gregate contained within a rigid cavity The effective 
mass M(ui) is complex valued; its real part reflects the 
inertial and clastic properties while its imaginary part re- 
flects the dissipative properties of the granular medium. 
We demonstrate how the features of M(u>) allow the 
study of some of the mechanisms of damping of acoustic 
modes, aging and elasticity in granular matter. 

Generally speaking, the real part of M{oj) exhibits a 
sharp resonance which we interpret in terms of an effec- 
tive sound speed. The imaginary part of M(w) shows a 
broad resonance peak which quantifies the attenuation 
of acoustic waves in the system. We observe significant 
changes in the stiffness and attenuation as a function of 
humidity. By monitoring the effective mass in time, we 
find a logarithmic aging effect in the resonance frequency 
as well as an increase of the damping upon humidifica- 
tion. These effects can be modeled as capillary conden- 
sation occurring between the asperities at the contact 
surfaces between the grains during humidity-dry cycles. 
We interpret this phenomenon in the context of a "trap 
model" of thermally activated liquid bridges. Our results 
suggest that, in the granular materials investigated in the 
present study, dissipation of acoustic energy is dominated 
by the asperities at the interparticle contact surfaces. In 
addition, humidity drastically affects the attenuation of 
the material through the capillary condensation of liquid 
bridges. 



Experiments. — A cylindrical cavity (of diameter 
D=25A mm) excavated in a rigid aluminum cup is filled 
with tungsten particles up to a height of 30.7 mm. We use 
tungsten particles due to their large density, which maxi- 
mizes the effects we arc studying. Each of these particles 
consists of four or five equal-axis particles, of nominal 
size 100 /xm, fused together (see inset of Fig. [I}. Such 
irregular arrangement produces low volume fractions (in 
this case 48%, in comparison with random close packing 
of spherical particles, 64%) and allows to study the loose 
packings of interest in the present work. To investigate 
the applicability of our results to other materials we per- 
form experiments with spherical glass beads as well as 
lead beads, see below. The cup is mounted on a B&K 
shaker and subjected to a vertical sinusoidal vibration at 
frequency lo = 2nf, with the frequency scan from sub 
kHz to 10 kHz. Since the bulk behavior depends on the 
contacts between the grains, the way of filling the gran- 
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FIG. 1: Effective mass of tungsten particles in a rigid cavity 
as a function of driving frequency. Inset: Scanning Electron 
Microscope image of the tungsten particles used in our exper- 
iments. The scale at the bottom of the image corresponds to 
300/im. 
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ular medium in the container is important. Prior to the 
measurement, we pour the grains into the cup and then 
shake the cup near its resonance frequency (see below) 
at large amplitudes up to 3g for 30 minutes to allow the 
compaction of the material. During the measurement, 
the granular medium is driven with cup acceleration of 
approximately 2 to 3 ms~ 2 . The amplitude of vibration 
is approximately 1 /im for the lower frequencies and de- 
creases to below 1 nm for the higher frequencies. 

Force F(oj) and acceleration a(oj) are independently 
measured with a force gauge and an accelerometer at- 
tached to the bottom of the cup. The signals are ana- 
logically filtered with a lock-in amplifier, which is syn- 
chronized with the driving signal. These filtered signals 
are then digitally recorded. The force gauge is sand- 
wiched between the shaker and the cup container. At 
each frequency, the force and acceleration are averaged 
over a 500 ms or longer time window. Prior to this time 
window, there is a 3 sec waiting time following the on- 
set of the driving frequency in order to allow the system 
to reach a steady state, and achieve accurate measure- 
ments. Since F(oj) and a(oj) are the net force exerted 
on and the response of the cup-granules assembly, their 
ratio, the effective mass, is independent of the shakers 
frequency response. Taking into account the mass of the 
empty cup, M c , we have M(oj) + M c = -j^y, where, 

M(oj) = M (oj) + iM (oj), reflects the partially in-phasc, 
M , partially out-of-phase, M , motion of the individ- 
ual grains, relative to the driving force. Thus M carries 
information about the inertia and elastic response while 
AI relates to attenuation of particle response. 

Fig. [T] shows the result for a frequency sweep of 
M(oj). We note that in the low frequency limit, M (oj — > 
0) — > 149g, the static mass of the grains. More interest- 
ingly, there is a relatively sharp resonance peak around 
/o = ojq/2w —1.8 kHz and a relatively broad tail that di- 
minishes with increasing frequency. M exhibits a broad 
peak centered at wq. We have conducted the experiments 
with other granular materials. Fig. [2Ja)-(b) shows the 
results for spherical lead beads of the size 180-250 /im 
and glass beads of 150 /mi. These results are qualita- 
tively similar to those of tungsten particles, in both the 
real and imaginary parts; the tungsten particles exhibit 
larger effects. 

We also applied the effective mass technique to liq- 
uid samples. Fig. [2}: shows a typical plot for the effec- 
tive mass of one such liquid (fluorocarbon fluid, density 
1800 kg/m 3 ) as measured by the effective mass technique. 
M (oj) shows a sharp resonance while the M (oj) curve is 
almost zero over the whole frequency range (except near 
resonance) indicating negligible dissipation. We can ex- 
tract a value for the effective compressional sound speed 
from the resonance frequency as v s (L) = 4£/o, where L 
is the depth of the medium. The values of density and 
sound speed for a variety of liquids, as measured with the 
effective-mass technique, are cross-plotted against those 
determined by more conventional means Q in Fig. 



The density of the liquids is measured using Archimedes 
principle and buoyancy (accuracy of 0.1%). The sound 
speed is measured with ultrasonic (centered at 500 kHz) 
through-transmission method (accuracy of 1%). The 
grain density of the tungsten sample is measured with 
a Micromeritics helium pyenometer (accuracy 1%) and 
the bulk density is calculated from the measured weight 
of tungsten grains filling a cup of a known volume. Vol- 
ume fraction is then estimated from the grain and bulk 
density. The results in Fig. (2J1 suggests that our tech- 
nique for measuring M(oj) is an accurate one. 

We investigate how the effective sound speed changes 
with the filling depth L, based on the peak frequency /o. 
The sound speed of the granular media is estimated based 
on the hypothesis of 1 /4 wavelength resonance frequency. 
Fig. [SJs shows the trend of greater speed with greater 
depth. The values we are reporting are of the same order 
of magnitude as those reported in the literature [l], [f| 
7]. In the range studied here, the linear relation v s (L) ~ 
L seems to hold up to L = 61 mm as seen in Fig. 

We may model the general features of Fig. Q] in 
terms of a simplified continuum model. The granular 
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FIG. 2: (a) Reai (M ) and (b) Imaginary part (M ) of ef- 
fective mass for spherical lead and glass beads, (c) Effective 
mass of fluorocarbon fluid, (d) Crossplot of measured values 
and book values of densities and sound speeds for various liq- 
uids: Methanol (o and □, first symbol density, second sound 
speed), Ethanol (+ and x), Fluorocarbon fluid (o and V), 
Silicon oil, (viscosity 10 mPa-s, A and *). (e) Effective speed 
of sound in tungsten granules as a function of filling depth, L, 
for constant cavity diameter, D — 25.4mm. (f) Continuum 
model prediction of M(u). 
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medium is considered to be a lossy fluid, with negligi- 
ble viscous effects at the walls. The acoustic pressure 
in the fluid is described by a one-dimensional model: 
P(z) = Asin(qz) + Bcos(qz). The boundary conditions 
are: (a) The displacement of the fluid at z = matches 
that of the cup, u(z — 0) = u CU p- (b) The acoustic pres- 
sure at the top of the fluid vanishes, P(z = L) = 0. Thus, 
P(z) = u cup paj 2 [siii(qz) — tan(gL) cos(qz)]/q. The force 
the cup exerts on the fluid is F = na 2 P(z = 0) and the 



acceleration of the cup is 
the effective mass: 



from which we obtain 



M(ui) = Tia 2 ^fpK /wtan(gL), 



(1) 



where K = Kq[1 — iuj£] is the (lossy) bulk modulus of the 
medium, and q = wy p/K is the complex wave vector 
with p the density of the medium. This model shows 
that resonance peaks, as seen in Fig. [2f, occur when L 
equals odd multiples of 1/4 wavelength of the acoustic 
wave in the medium, with the main resonance peak at 
/o = f / (4L) yj Kq I p. Encouraged by these results, we 
interpret the main resonance in Fig. [1] as being a 1/4 
wavelength resonance of the compressional sound speed. 

Next, we use the features in M(u>) to investigate the 
dissipation in the granular medium. The dissipation 
mechanisms can be grouped into two categories: 1) Local 
damping throughout the contact network that includes 
asperity deformation and bulk attenuation Q, and wet- 
ting dynamics in the liquid bridges caused by relative 
motion of neighboring particles [9(. Wetting dynamics 
also results in greater contact stiffness [l(| • If large inter- 
particle tangential displacements occur then, there exist 
dissipation through Coulomb friction as well. 2) Global 
damping due to drag caused by particle motion in a sur- 
rounding viscous medium [Til ] . 

Next, we perform further experimentation by changing 
the ambient conditions. Ambient humidity can modify 
the dissipation and the stiffness at the grain to grain sur- 
face contacts [Toj . For instance, small amounts of liquid 
greatly affects the maximum angle of stability in sand- 
piles [12| , while capillary action of the interstitial li quid 
results in increased cohesion [l(| , modified friction [13| 
and aging in the angle of repose of granular materials 
[l4| . We place the effective mass apparatus in a chamber 
where humidity is controlled by first using an aqueous 
saturated salt solution [l5| and subsequently a dessicant. 
We measure the relative humidity as the ratio of the par- 
tial pressure of water to the saturating vapor pressure 

Hk = Pvap/Psat- 

We follow a humid-dry cycling process as follows: we 
start from room humidity = 26% at t < 0, and then 
we humidify at t = to -Hr = 75% until t = 126 hr, and 
then we dry to H-g, = 0% until t = 174 hr. Fig. [3] shows 
the effective mass measurements under the humidity con- 
trolled environment carried out at various time intervals. 
We observe a very slow shift in resonant frequency during 
the humidification cycle, followed by a further positive 
shift upon drying. The effect is more important when we 
repeat the same cycling but to 95% humidity. 
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FIG. 3: Imaginary part of M(f) measured under humidity 
controlled conditions (the real part shows similar behavior). 
The solid lines represent measurements under Hb, = 75% and 
the dashed lines are for the dry state, Hb. = 0%. Top in- 
set: Frequency shift vs time for _Hr=75% (o) and 95% (□). 
Bottom inset: Damping rate (r) vs time for Hn—75%. The 
dashed lines in the insets are logarithmic fits. The up and 
down arrows indicates the large variations in stiffness and 
damping upon drying, respectively. 



During humidification, we observe an aging phe- 
nomenon as evidenced by the approximate logarithmic 
behavior of the shift in resonant frequency Afo(t) (top 
inset of Fig. 0), A/ (t) - C(H R ) log(i/t<). where C(H R ) 
is a constant that increases with humidity, C(75%) = 292 
Hz and C(95%) = 494 Hz and U « 2 hr for both mea- 
sures. 

Furthermore, the humid system shows a large increase 
in the attenuation as illustrated by the increase in damp- 
ing rate V (bottom inset of Fig. [3]) . Based on the exper- 
imental effective mass data, we search for poles in M(Q) 
for Co in the complex plane. Once the complex resonance 
frequency (uiq = lu + iu> ) is found, the damping rate is 
given by T = —lu q . The damping rate at ambient humid- 
ity is 300 s _1 , which rises to 1200 s _1 under prolonged 
humidification (see bottom inset of Fig. [3]). The damping 
data suggests logarithmic increase with time. 

A possible interpretation of the observed behavior sug- 
gests that contact forces through capillary liquid bridges 
arc important for damping, stiffness and aging. As 
demonstrated in [l4| . the humidification process leads 
to the formation of capillary bridges around the rough- 
ness asperities at the interparticle contact surfaces. The 
aging effect observed in the logarithmic variation of the 
frequency shift with time can be interpreted as originat- 
ing from capillary condensation. The key quantity to 
calculate is the time dependence of the number of wetted 
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asperities n a {t). Capillary condensation at the surface 
asperities can be seen as a thermally activated process 
following a model proposed in The energy barrier 

to form a liquid bridge of cross sectional area A between 
two asperities separated by a gap h is AE = AppiAh, 
where pi is the molecular density of the condensing liq- 
uid, Ayit = U-qT ln(l I Hr) is the undersaturation of the 
chemical potential, and T the room temperature. The 
activation process follows an Arrhenius dynamics such 
that the time to create a bridge is t = t exp(AE/k B T), 
with t the microscopic condensation time. 

Due to the roughness of the particles we expect a dis- 
tribution of gaps P(h), which leads to a distribution of 
energy barriers, P(AE). The calculation of the probabil- 
ity distribution P(t) is then analogous to obtaining the 
probability of the time to escape from a trap with energy 
barrier AE for a thermally activated particle in a random 
energy potential. Such a process corresponds to the "trap 
model" used to describe the phenomenology of glasses 
(Til. The probability P(t) satisfies P(t)cIt = P(h)dh. 
While the probability P(h) is unknown, we expect it 
to be a bounded distribution determined by the surface 
roughness of the particles. The two most common cases 
lead to the same result: for an exponential or Gaussian 
distribution of gaps, wc find P(t) ~ r _1 . 

The number of activated wetted asperities at time t 
is n a (t) oc f t P(r)dr. Then n a (t) w log(t/t ). Liquid 
bridges induce a significant cohesion between the parti- 
cles which for loose grains under gravity and large hu- 
midity outweigh the weight of the grains. Therefore, the 
logarithmic increase in n a leads to an additional stiffness 
with the concomitant increase in the compressional sound 
speed. Since the resonance frequency is proportional to 
the stiffness of the contact, the logarithmic aging can be 
interpreted as arising from the wetted asperities. 

Upon drying the particles to £/r, = 0% after humidifi- 
cation, the physics is very different from the humid pro- 
cess. As shown by the experiments of fioj . the capillary 



rings dry out and the asperities are sucked down by the 
increasing capillary pressure in the evaporating film lead- 
ing to an even larger increase in stiffness of the contacts. 
Thus, in the dried out state there is much more solid 
on solid contact and the sound speed is increased. This 
increase in stiffness upon drying results in even larger 
positive shifts in resonant frequency as seen in the top 
inset of Fig. [3] 

We also find important changes in the dissipation in 
the dry state. Upon drying, the decrease in overall atten- 
uation is evident from the much sharper resonance in M 
(we find a very large mass at resonance M (/o) ~ 3500 
g). Accordingly, the damping rate decreases to Y — 50 
s" 1 in the dry case from Y = 1200 s" 1 at H R = 75% 
(see bottom inset of Fig. [3]). In the driest state, no liq- 
uids bridges are left and the dissipation is purely through 
asperities deformation with the concomitant lower damp- 
ing strength. These results suggest that, for the granular 
media studied here when humid conditions prevail (in- 
cluding ambient), liquid bridges provide larger attenua- 
tion of acoustic waves than the viscoclastic dissipation in 
the bulk of the particles or viscous global damping. 

In conclusion, we have demonstrated the usefulness of 
the effective mass technique to study the dissipation of 
acoustic waves, stiffness and aging in granular media. Ex- 
periments and theoretical models suggest that capillary 
condensation at the interparticle contacts is an impor- 
tant mechanism for dissipation and aging. The kinetics 
of aging can be modeled in terms of a trap model for the 
formation of liquid bridges at the asperities. The large 
variations in dissipation found in this study demonstrate 
the conditions for effective particle damping and are rel- 
evant to a variety of applications for optimizing attenu- 
ation of structure-borne acoustic waves. 
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